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Abstract 

The kj_ dependent gluon distribution for protons, pions and photons is calculated numer- 
ically accounting for the resummation of small x effects due to the Lipatov equation. It 
is represented by a convolution of a gluon density in the collinear limit and a universal 
function Q(x,k 2 ,fi) for which an analytic expression is derived. 
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1 Introduction 



In the small x range new dynamical effects are expected to determine the behaviour of structure 
functions. The evolution of parton densities is effected by terms due to non strong k± ordering 
and eventually by screening contributions. This behaviour is expected to determine not only the 
structure functions of the nucleon but also the parton distributions of the pion and the photon. 
A description of contributions of this type requires to generalize the factorization of the hadronic 
matrix elements into coefficient functions and parton densities in which the k± dependence is 
not integrated out [1]. 

This factorization covers the case of collinear factorization in the limit that the k 2 depen- 
dence of the coefficient function is neglected. The k± dependent gluon density accounts for the 
resummation of small x effects. In the present paper we will consider those due to the Lipatov 
equation only. Since this equation behaves infrared finite no other singularities will emerge than 
in the case of mass factorization. The collinear singularities are delt with in the same way in the 
case of k± factorization. 

In the present paper the k± dependent gluon distribution is calculated for the case of the 
scheme [2, 3] for the proton, pion, and the photon. The k± dependent distribution k 2 } ji) can 
be represented as the convolution of the gluon density in the collinear limit g(x } ji) and a function 
Q(x } k 2 } fi) for which an analytic expression will be derived. A numerical comparison of the k± 
dependent gluon densities accounting for the solution of the Lipatov equation and in the double 
logarithmic approximation (DLA) is given for the different hadrons and the photon. It may be 
interesting to try to unfold these distributions using different observables (F 2 , Fl, cr 2 Je t, ^j/tp) 
and to compare the behaviour of this distribution for different hadrons (p, ir) and the photon to 
obtain a conclusive understanding of this quantity. Here we focus on the quantitative behaviour 
of k 2 } ji) as a function of x and k 2 to investigate first the size of the novel effects in comparison 
with standard ones. 



2 k±_ Factorization and the k±_ dependent gluon distri- 
bution 

The factorization relation for an observable Oi(x } fi) reads 

Oi(x } fi) = J dk 2 (70i(x, k 2 } ji) ® k 2 } ji) (1) 

where aoi(x, k 2 } ji) and $(x,A; 2 ,/i) denote the k 2 dependent coefficient function and parton 
density 1 , respectively. Eq. (1) can be rewritten as [2, 3] 

Oi(x, fi) = (Jq.(x, fi) (x) G(x, fi) + / dk 2 ^aoi(x, k 2 , /i) — (Jq.(x, k 2 , /i) (2) 

J 

with <Tq = lim fc 2_j. aoi(x, k 2 } ji). The first addend in (2) describes the conventional con- 

tribution due to collinear factorization. The second term contains the new contributions. Note 
that k 2 } fi) starts with terms oc a s . It has therefore not the interpretation of a probability 
density and may even become negative. 

1 We will consider the gluon density in the present paper only. 



As shown in [2] the k± dependent gluon distribution associated to eq. (2) reads in moment 
space 

1 /p\ 7c ( J '~ s ) 

k 2 ,p) = 7 c (j, a s )— f-J g(j, p) (3) 

where p denotes a factorization scale, a s = N c a s (p)/ir , and g(x } p) is the gluon density. Eq. (3) 
accounts for the small x behaviour due to the Lipatov equation. Here 7 c (j, ct s ) 1S the solution 
of the eigenvalue equation of the homogeneous Lipatov equation 2 

= x(7cU, «*)), x(7) = 2V>(l)-V>(7)-V>(l-7)- (4) 



In x space the k± dependent distribution is given by the convolution 

k 2 , p) = Q(x, k 2 , p) (g) g(x, p), (5) 

correspondingly, with 

\ dk 2 §(x,k 2 ,p) = 5(1 - x). (6) 

JO 

The function Q(x } k 2 } p) is universal and can be calculated numerically by a contour integral 
in the complex plane over the first factor in eq. (3). Since the solution of eq. (4) is multivalued 
the Mellin inversion to x space requires to select the branch in which for asymptotic values of 
j e C 7 C approaches the perturbative result 7 c (j, ct s ) ~ ®s/(j — 1) for small values of a s . 

We solved eq. (4) under this condition numerically using an adaptive Newton algorithm. The 
solution is characterized by three branch points. Their position in the p plane may be determined 
using the solutions of (7) in eq. (4) 

l sin [ttz) 

Obviously z = 7* = 1/2 is one of the solutions of (7) leading to p\ = 4 In 2. The two other 
solutions are 

7c s (2 3) = -0.4252 ± 0.4739 i 

p s (2j3) = -1.4105 ± 1.9721 i. (8) 

This agrees with a result in ref. [5]. As shown in figure 1 the smooth dependence of "f(p) 
for Rep > 4 In 2 changes if Rep < Rep\. I1T17 becomes discontinuous at Imp = and Re7 
developes a 'ridge' with Re7 = 0.5 at Imp = 0. In figure 2 and 3 the real and imaginary 
part of 7 C are illustrated for smaller values of Rep. For the real part the 'ridge' persists until 
Rep ~ —1.3 and turns into a flat form with Re(7 c ) ~ —0.5 for | Imp | < 1.5. At the same time 
the imaginary part of 7 C becomes continuous again. We have also shown the contour of Re7 c 
and Im7c for Rep = Rep^>3)- The branch points are clearly seen as edges in Re7 c and tips in 
Im7c. The curvature of Re, lm(7 c ) changes its sign upon passing the branch point at Rep = const, 
discontinuously. At slightly larger values of Rep this change is continuous still (cf. figure 2,3). 



2 In several recent approaches [4] to describe k±_ dependent gluon distributions phenomenological Ansatze were 
used based on solutions of inhomogeneous Lipatov equations. Note that these descriptions are not related to 
eq. (2) and [2, 3]. 
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Figure 4: The kj_ dependent gluon distribution of 
the proton as a function of k 2 and x. Full lines: 
complete solution; dashed lines: solution in DLA. 
For the input distribution g(x,fi) the parametriza- 
tion [8] (LO) was used. 



3 An analytical solution for ^(x,/c 2 ,/i) 



The integration contour for the Mellin transformation of A; 2 , p) to A; 2 , p) for j eC has to 
be situated outside the range of the singularities of 7 C . One may expand 7 c (j, a s ) into a Laurent 
series over p 

oo 

7c(j, «*) = Y,9W~ l (9) 

in this range. 

The coefficients g\ are given in [6] up to / = 20 in analytical form extending an earlier result [7]. 
For small values of | Imp | (|lm/?| < 2) the truncated Laurent series leads to an oscillatory behaviour 
and eq. (9) does no longer serve to be an appropriate description of 7 C (cf. [6]). 

Using (9) a corresponding expansion may be performed for 

k 2 g(j, k 2 ,p) = 7 c (j, a s ) exp[7 c (j, a s )L] (10) 

with L = \n(k 2 / p 2 ). For the single terms of the Laurent series in p the Mellin transform can be 
carried out analytically. 

Here it is important to expand the exponential in eq. (10) in such a way that the lowest order 
term in a s of 7 C is kept in exponential form. One obtains 

/— r777T7A , j fr ^{^log(l/x)\ {v - 1)/2 



k 2 Q(x,k 2 ,p) = ^I (2Ja s \og(l/x)L) +^STd v (L) 



L 



x ^a s log(l/x)^ ,L>0. (11) 

The coefficients d^L) are given in ref. [6]. Up to v = 20 they contain at most terms oc L 4 . The 
first term in eq. (11) denotes the Green's function in DLA. 
For L — > (11) takes the form 



S('-i) 

and for L < (i.e. A; 2 < p 2 ) one has 



n '-1 

(12) 



k 2 g(x,k 2 ,>L) = ^J (2JaJog(l/x)\L\)+^f:d4L)^ OslOg{1/x) ^ ^ 



x J^-i ^a s log(l/x)|L|) . (13) 
Thus for A; 2 — > damped, oscillating modes are obtained which vanish faster than l/l^l" 1 / 4 . 



4 Numerical Results 

The k± dependent gluon distribution $(x, A; 2 , p) (scaled by A; 2 ) is shown in figure 4 as a function 
of x and A; 2 for p 2 = 20 GeV 2 for the proton referring to the parametrization of ref. [8] to describe 
g(x,p). The complete solution eq. (5) is larger than the DLA result for A; 2 ~ p 2 at x ~ 10~ 3 by 
10 to 15% while for A; 2 — > smaller values are obtained. 
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Figure 5: The kj_ dependent gluon distribution Figure 6: The kj_ dependent gluon distribution of 

of the pion as a function of k 2 and x. Full the photon scaled by 1 / clqed as a function of 

lines: complete solution; dashed lines: solution k 2 and x. Full lines: complete solution; dashed 

in DLA. For the input distribution g(x,fi) the lines: solution in DLA. For the input distribution 

parametrization [9] (LO) was used. g(x,fi) the parametrization [10] (LO) was used. 

At larger values of x the complete solution approaches the DLA result. For k 2 — > k 2 } fi) 
vanishes. Since the DLA result is proportional to J (2^a s log(l/x) log \k 2 /fi 2 \) for k 2 — > a 
damped oscillatory behaviour is obtained in this approximation. The complete solution, on the 
other hand, behaves monotonous in the whole kinematical range. 

The calculation of Q(x } k 2 } fi) in eq. (5) by numerical Mellin inversion using the numerical 
solution of eq. (4) is rather time consuming compared to the convolution of the analytical solution 
(sect. 3) with the input distribution g(x } fi). We compared both methods and found that the 
representation given in the previous section leads to a relative error of less than 0.002 using an 
expansion up to O(a 2 s ). 

The same calculation has been performed for the k± dependent gluon distribution of the 
pion and the photon using the distributions [9, 10] (LO) as input parametrizations. The ratio of 
the complete solution and the solution in DLA is of similar size than in the case of the proton. 
Comparing the functions k 2 , /i) and (Xq 1 ed $ 1 (x, k 2 , /i) the latter one turns out to be larger 

at virtualities k 2 ~ jj 2 . 

Since shape and size of the complete solution and in DLA are rather similar very precise 
measurements are required to establish the non-DLA contributions at small x both in the case 
of hadrons (p, ir) and photons. 



5 Conclusions 

We have calculated the k± dependent gluon density for the proton, pion, and photon numerically 
in leading order using the BFKL equation. A consistent treatment of observables is possible in 
the scheme [2, 3]. 



The Green's function Q(x } k 2 } fi) was found both numerically and by a perturbative analytic 
expression expanding the complete solution up to 0(a 2 s °). Both representations agree better 
than 0.002 after convoluting with the respective input distributions g(x } fi). 

The effect of the non-DLA terms in k 2 } fi) is of (9(10. ..15%). To reveal these contribu- 
tions requires very accurate measurements in the small x range in the case of every observable 
being sensitive to the gluon density widely independent of the choice of the target. 
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represented by a convolution of a gluon density in the collinear limit and a universal 
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